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Abstract 

Was the spread of agropastoralism from the Fertile Crescent throughout Europe influenced by extreme climate events, or 
rsq was it independent of climate? We here generate idealized climate events using palaeoclimate records. In a mathematical 
y—i model of regional sociocultural development, these events disturb the subsistence base of simulated forager and farmer 
^^ societies. We evaluate the regional simulated transition timings and durations against a published large set of radiocarbon 
^^ dates for western Eurasia; the model is able to realistically hindcast much of the inhomogeneous space-time evolution 
bJOof regional Neolithic transitions. Our study shows that the consideration of climate events improves the simulation 
^ of typical lags between cultural complexes, but that the overall difference to a model without climate events is not 
"^>, significant. Climate events may not have been as important for early sociocultural dynamics as endogenous factors. 

' Key words: Europe, climate events, extreme events. Neolithic transition, adaptation, modeling 
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1. Introduction 

Between 10 000 and SOOOcalBC, western Eurasia saw 
enormous cultural, technological, and sociopolitical changes 
with the emergence of agropastoralism, permanent settle- 
ments, and state formation (Barker, 2006). Human pop- 
ulation experienced a dramatic increase (Bocquet-Appel, 
2008; Gignoux et al., 2011), and people, plants and anim- 
als moved or were moved great distances (e.g., Zohary and 
Hopf, 1993). 

While the Holocene possibly defines the start of ma- 
jor anthropogenic global environmental change (Lemmen, 
2010; Kaplan et al., 2011), it also marks the period where 
climatic shifts could have affected human subsistence more 
severely than ever before: reduced mobility after invest- 
ments in settlement infrastructure most likely increased 
the sensitivity of the novel farmers to environmental alter- 
ations (Janssen and Scheffer, 2004). There remains, how- 
ever, considerable uncertainty on whether and how climate 
instabilities had infiuenced the development and spread of 
agropastoralism in Eurasia (Berglund, 2003; Coombcs and 
Barber, 2005). 

1.1. Origin and spread of western Eurasian farming 

The Neolithic originated most probably in the Fer- 
tile Crescent, between the Levantine coast and the Zagros 
ridge. In this region, almost all European food crops and 
animals — wheat, barley, cattle, sheep, pigs — had been do- 
mesticated and inserted into a broad spectrum of foraging 
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practices during the tenth millennium calBC (Flannery, 
1973; Zeder, 2008). Neohthic (farming based) hfe style 
emerged not before the 9th millennium BC in this core re- 
gion (Rosen and Rivera-CoUazo, 2012), and expanded to 
Cyprus by 8500calBC (Peltenburg et al., 2000); around 
7000calBC, agropastoralism appeared on the Balkan and 
in Greece (Perles, 2001). Propagating in a generally north- 
western direction, agropastoralism finally arrived after 4000 
cal BC on the British isles and throughout northern Europe 
(Sheridan, 2007); in a western direction, the expansion 
proceeded fast along the Mediterranean coast to reach the 
Iberian peninsula at 5600 cal BC (Zapata et al. , 2004) . 

1.2. Transitions and climate 

It has been argued that a precondition of agriculture 
was the relatively stable environment of the Holocene (Feyn- 
man and Ruzmaikin, 2007), and that only in this stable 
environment active cultivation and establishment of infra- 
structure such as fields and villages was favored (van der 
Leeuw, 2008). Within its relative stability, however, the 
Holocene climate exhibited variability on many spatial and 
temporal scales with pronounced multi-centennial and mil- 
lennial cycles (Mayewski et al., 2004; Wanner et al., 2008). 
In addition, non-cyclic anomalies have been identified (e.g., 
Wirtz et al., 2010), most prominently the so-called 8.2 
and 4.2 events (around 6200 and 2200 cal BC, respectively, 
von Grafenstein et al. 1998; Cullen et al. 2000). Although 
the regional scale and intensity of the 4.2 event has been 
strongly questioned (e.g., Finne et al., 2011), the event had 
evoked the formulation of hypotheses on the connection 
between climatic disruptions and societal collapse (Weiss 
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et al., 1993; DeMenocal, 2001). Similarly, the globally doc- 
umented 8.2 event has been linked to the abandonment 
of many settlements in the Near East and simultaneous 
appearance of new village structures in southeast Europe 
(Weninger et al., 2005). 

It might be coincidental that the 8.2 and 4.2 events 
define the time window of the Neolithic expansion in Europe, 
but the general view that environmental pressure on early 
Neolithic populations may have stimulated outmigration 
has been put forward since long (Childe, 1942). Dolukhanov 
(1973), Gronenborn (2009, 2010), or Weninger et al. (2009) 
suggest that climate-induced crises may have forced early 
farming communities to fission and move in order to escape 
conflicts. Berger and Guilaine (2009), to the contrary, see 
the role of climate events rather in creating opportunities: 
the rapid farming expansion into the Balkan could have 
been stimulated by an increase of natural fires after the 
8.2 event, which opened up the formerly forested land- 
scape. 

1.3. How sensitive was the Neolithization to climate? 

The relevance of climate variability and external trig- 
gers for prehistoric agricultural dynamics has been severely 
questioned (e.g., Erickson, 1999; Coombes and Barber, 
2005). Alternative theories of the Neolithic transition un- 
derline the agency of early societies (Shanks and Tilley, 
1987; Whittle and Cummings, 2007). On the other hand, 
the development of technological, social, and cultural com- 
plexes can hardly be thought to evolve independently of 
their variable environments; and the spatio-temporal im- 
print of the Neolithization in Eurasia requires a geographic 
approach which resolves how people and/or goods and 
practices migrated over long distances. Berglund (2003), 
e.g., suggested a stepwise interaction between agriculture 
and climate but found no strong links for northwest Europe. 

The dispersal of agriculture into Europe has long been 
mathematically formulated based on Childe's (1925) ob- 
servation on the spatiotcmporal distribution gradient of 
ceramics that Ammerman and Cavalli-Sforza (1971) for- 
mulated as the 'wave of advance' model. This simple — and 
also the later more advanced ones (Ackland et al., 2007; 
Galeta et al., 2011; Davison et al., 2006) — diffusion models 
received support from linguistic (e.g. Renfrew, 1987) and 
archaeogenetic work (e.g. Balaresque et al., 2010). The 
dispersal of agriculture in these models occurs concentric- 
ally, and can be modulated by topography and geography. 
This dispersal model is not able to describe the inhomo- 
geneous spatiotcmporal distribution of radiocarbon dates, 
which are, e.g., apparent in regionally different stagnation 
periods ('hypothese arythmique', Guilaine, 2003; Rasse, 
2008; Schier, 2009). 

Stagnations are visible in the simulation by Lcmmcn 
et al. (2011), who integrate endogenous regional sociocul- 
tural dynamics with the dispersal of agriculture. Their 
approach connects social dynamics — as optimally evolving 
agents — to regionally and temporally changing environ- 
ments; in addition, they account for the spatio-temporal 



spread of populations and technological traits. Their Global 
Land Use and technological Evolution Simulator (GLUES) 
has proven to produce realistic hindcasts of the origin and 
distribution of agropastoralism and concomitant cultures 
around the globe (Wirtz and Lemmen, 2003; de Vries et al., 
2002), for Eastern North America (Lemmen, 2012), the In- 
dus valley (Lemmen and Khan, 2012), and western Eurasia 
(Lemmen et al., 2011). Using GLUES and a globally syn- 
chronous climate forcing signal, Wirtz and Lemmen (2003) 
found a general delay of the simulated regional Neolithic 
due to climate fluctuations; at a global scale, differences in 
hindcasted socio-cultural trajectories proved to be largely 
independent of temporal disruptions. 

We here use temporal disruptions that are defined as 
excursions of a climate variable far from the local mean 
climate, i.e. extreme climate events; we do not consider 
rapid climate shifts that abruptly alter the climate mean 
state (e.g., Dakos and Scheffer, 2008). The hypothesis that 
extreme climate events had significant impacts on the Neo- 
lithization of Europe is critically examined: we employ 
GLUES as a deductive tool to reconstruct the Neolithic 
transition in Europe and evaluate the simulated recon- 
struction against the radiocarbon record of Neolithic sites 
in two experiments: (1) one including climate events, rep- 
resented by a pseudo-realistic spatially resolved climate 
event history for the period 9500-3000 calBC; and (2) an- 
other without climate events. 

2. Material and Methods 

2.1. Reconstructing climate event history 

We used a data collection of 134 globally distributed, 
high- resolution (< 200 a) and long-term (> 4000 a) pa- 
laeoclimate time series collected from public archives and 
published literature. The collection only contains studies 
where the respective authors indicated a direct relation 
to climate variables such a precipitation, temperature, or 
wind regime (e.g. Bond, 1997; Wick et al., 2003; Chapman 
and Shackleton, 2000; Gasse, 2000). A large part of this 
data set (122 time series) was previously analyzed by Wirtz 
et al. (2010) for extreme events; a complete overview of 
time series in this collection is provided in the supporting 
online material (Table SI). Due to the different types of 
proxies originating from both marine and terrestrial sites 
(mostly (5"'^*0, see Table 1) the relation to climate vari- 
ables is often ambiguous, also in sign. This ambiguity 
does not affect our analysis, as we are only interested in 
the spatio-temporal characterization of extreme events: a 
drastic excursion from a climate mean state stressed re- 
gional habitats and human populations regardless of its 
direction. 

Our data set comprises 134 palaeoclimate time series, 
all of them long-term and high-resolution, and provides 
the best spatial and temporal coverage of any study we are 
aware of. Previous collections used 18, 50, 60 or 80 records 
(Wanner et al. 2008; Mayewski et al. 2004; Holmgren et al. 



2003; Finne et al. 2011, respectively), mostly limited to 
the last 6000 years. The coverage we use here is sufficient 
to represent climate variability in almost all land areas 
of the world (with sparsest regional coverage in central 
Australia, Saharan Africa, and Northern-Central Eurasia) , 
considering the spatial coherence of climate signals within 
1500 km distance found by Wirtz et al. (2010) for their 
similar data set. 

From the global dataset, 26 time series are located in 
or near our focus area western Eurasia (Table 1). For these 
time series we analyzed the non-cyclic event frequency ac- 
cording to the procedure in Wirtz et al. (2010): time series 
were dctrended with a moving window of 2000 years and 
smoothed with a moving window of 50 years, then normal- 
ized (Figure lb). Events were detected whenever a time 
series signal exceeded a confidence interval with threshold 
p — 1 — 1/n, where n is the number of data points (Thom- 
son, 1990), and where each event is preceded or followed 
by a sign change in the time series. 

For each simulation region, events from spatially over- 
lapping or nearby proxy locations were used to construct 
an aggregated event time series specific to this region (Fig- 
ure la-c): (1) a Gaussian filter with r = 175 a (corres- 
ponding to the dating uncertainty in many records) stand- 
ard deviation was applied to each event; (2) the distance 
of the proxy location to the simulation region was used 
to assign exponentially decreasing weights to each event; 
(3) all event densities were summed and filtered with a 
running mean with window t; (4) all events representing 
single anomalies above the mean event density were used 
for further analysis, irrespective of their magnitude. On 
average, seven events were detected per time series. 

From this analysis, we obtain a pseudo-realistic data 
base of spatially and temporally resolved climate events. 
By using this idealized approach, we show a way to over- 
come issues raised by, e.g., Schulting (2010) on chrono- 
logical resolution and spatial representation problems of 
individual records. In our study, we consider the impact 
of the number of extreme events and their spatial pat- 
terning rather than single chronologies. This idealized — 
or potential — climate events database thus allows us to 
go forward with analyzing the climate-human relationship 
while the reliability of individual palaeoclimate reconstruc- 
tions is still being questioned. 

2.2. Global Land Use and technological Evolution Simu- 
lator 

The Global Land Use and technological Evolution Sim- 
ulator (GLUES, Wirtz and Lemmen 2003; Lemmen 2010; 
Lemmen et al. 2011) was developed to study how dif- 
ferences in cultural trajectories at the region scale can 
be attributed to the specific adaptation of local societies. 
GLUES mathematically resolves the dynamics of local hu- 
man populations' density and characteristic sociocultural 
traits in the context of a changing biogeographical environ- 
ment. One of the characteristic traits represents available 
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Figure 1: Generation of idealized climate events, (a, top) Map 
with example region and nearby locations (colored triangles) where 
long-term and high-resolution palaeoclimate information is available; 
concentric circles indicate the distance from the simulation region. 
(b) Associated palaeoclimate time series, detrended and normalized; 
for each time series, the confidence level p is indicated and visual- 
ized as the width of the background shading. Climate extremes are 
identified outside this confidence interval and are highlighted, (c, 
bottom) Contribution of individual extreme events from palaeocli- 
mate time series (weighted by distance and color-coded as above) 
to the probability distribution of event occurrence in the simulation 
region. The dotted line shows the number of time series contributing 
to the event generation for this region. 



technologies, a second one the share of farming and herd- 
ing activities, and the third one the number of established 
agropastoral economies. Trait characteristics adapt ac- 
cording to a growth-benefit gradient dynamics (e.g. Smith 
et al., 2011, for ecological applications). Traits are further 
exchanged between simulation regions by information dis- 



Table 1: Table of 26 palacoclimate time series and associated timing of extreme climate events. This is a subset (relevant for the Western 
Eurasian focus area) of the entire global collection used to generated extreme climate events; the full dataset is shown in the supplementary ma- 
terial Table SI. Abbreviations: Lk=Lake, Cv=Cave, SST=sea surface temperature, T=temperature, T7=July temperature, P=precipitation, 
GSD=grayscale density, Ti=Titanium content, and HSG=hematite stained glass. 
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persal or migration. The model is described in more detail 
in the supporting online material. 

Our climate event data base is used to impose sudden 
reductions (/) of the land utility whenever a climate event 
occurs in a simulation region. If a local simulated popula- 
tion experiences reduced food yield during a climate event, 
people attempt to mitigate the crisis by moving to adjacent 
regions. In this respect, climate events act as a stimulus 
to migration in the model. If outmigration is not pos- 
sible, population density declines; with severe population 
decline, also some of the technologies are lost (Lcmmen 
and Wirtz, 2010). 

The simulation is started at 9500simBC^. All of the 
685 biogeographically defined regions are initially set with 
farming activity at 4% and established agropastoral com- 
munities at 25%, what represents a low density Mesolithic 
technology population and a broad spectrum foraging life- 
style with low unintentional farming activity. We define a 
local population as Neolithic when the share of agropastor- 
alists is larger than the share of foragers — regardless of its 
technology, economic diversity, or population density (Fig- 
ure 2). Experiments are performed with different impact 
strength (/) of extreme events on land utility to assess cli- 
mate related sensitivities; here we compare the simulation 
without climate events (/ = 0%) to simulations where the 
utility reduction was / = 10% to 100%; results are shown 
for / — 40%, a value that represents a moderate but not 
excessive impact. 

2.S. Reference data 

Our reference is a sub-set of the comprehensive data 
collection by Pinhasi et al. (2005), who used site data 
provided by the United Kingdom Archaeology Data Ser- 
vice, the Central Anatolian Neolithic e- Workshop (CANeW), 
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Figure 2: Simulated climate events and evolution of farmer share 
for an example region, (a, top) Declining resource availability at 
times where the palaeo climate record indicates climate extremes, (b, 
bottom) Transition from hunting-gathering to agropastoralism in the 
simulation without (grey) and with (black) climate events. Both the 
timing of the transition as well as the duration of the transition 
(defined as the time needed in the simulation to increase the farmer 
share from 4.1% to 95%) are prolonged by the introduction of climate 
events in this example region. 



the radiocarbon CONTEXT database, and the Radiokar- 
bondaten Online (RADON) database. In their compila- 
tion, Pinhasi et al. included only sites with small dating 
uncertainty (< 200 a) and report dates based on calibra- 
tion of original ^''C measurements with CalPal 2004. The 
data set contains dates distributed across Western Eurasia 
and was recommended by its authors for comparison with 
simulation studies (see Figure 3 for dates and locations). 
For numerical comparison with the simulation, we calcu- 
late the timing of the Neolithic onset from the radiocarbon 
date of all sites located within each simulation region. As 
there is no standard procedure for calculating such an area- 



averaged onset (Miillcr, 2000), we us the average date. 

3. Results 

In the reference GLUES simulation including climate 
variability, farming originates in the Fertile Crescent and 
southeast Europe in the 7th millennium simBC and pen- 
etrates into Europe in a northwest direction (Figure 3). 
By 3500simBC all of continental Europe has converted to 
farming as the predominant subsistence style. 

3.1. Expansion of agropastoralism 

The initial development progresses slowly and at a low 
level. It begins during the 67th century simBC in the 
Levant and Greece^, followed by the central Balkan (66th 
century). From these centers, farming spreads to the west- 
ern Black Sea coast during the 62nd century, and is present 
throughout the entire Balkan, Anatolia, and Mesopotamia 
by the 59th century. In these areas, the northward spread 
of farming stagnates, while farming activity intensifies up 
to the 56th century on the Balkan, and while farming ex- 
pands to Italy and the eastern Black Sea coast. 

From the 55th to the 59th century simBC, a rapid 
expansion through central Europe — or what is archaeolo- 
gically seen as the Linearbandkeramik area — is simulated, 
such that farming is the major subsistence style through- 
out central and southeastern Europe, including the the 
northern Black Sea coast. From the 46th century, farm- 
ing emerges along the Baltic Sea coast, and has become the 
dominant subsistence style throughout the North European 
plains, Poland, and Denmark by the 42nd century. In 
France and England, first farming is evident from the 39th 
century. 

Independently, farming also originates in North Africa 
at the Gibraltar Strait in the 61st century simBC, from 
where it spreads into Morocco and the Iberian peninsula, 
penetrating southern Spain and Morocco by the 54th cen- 
tury. The expansion into northern Iberia stagnates until 
the 49th century; and slowly connects along the Medi- 
terranean coast to the Fertile Crescent expansion branch 
by the 46th century. All of of the Iberian peninsula has 
converted to farming by the 42nd century. By the 34th 
century, all of continental western Eurasia and England 
rely almost exclusively on farming. An animation show- 
ing the expansion of agriculture into Europe is shown as a 
supplementary Movie SI. 

While the simulation captures the eastern route of the 
Neolithic into Europe (via Greece, the Balkan, Hungary, 
then north- and westward, see Rasse 2008 for an over- 
view of routes) quite well, it fails to simulate the Medi- 
terranean route (Cyprus, Sicily, French coast, then north- 
ward); this was attributed by Lemmen et al. (2011) to 
the lack of sea transport in the current model. A western 
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route into Europe is suggested by the model emanating 
from the Strait of Gibraltar; apart from some zooarchao- 
logical evidence (Andcrung et al., 2005), this route has not 
been confirmed by archaeology (Gronenborn, 2009; Rasse, 
2008). 

The regional timing of agropastoralism is contrasted 
with the median radiocarbon dates of Neolithic sites com- 
piled by Pinhasi et al. (2005) (Figure 3). From this synop- 
tic, time-integrated perspective, the simulated centers of 
agropastoralism in the Fertile Crescent, in northern Greece 
and at the Strait of Gibraltar are evident, as well as the 
southeast to northwest temporal gradient of the Neolithic 
transition. The general pattern of the simulated transition 
resembles the pattern that can be seen from the radiocar- 
bon dates; locally, many dates deviate from the large-scale 
simulation or disagree with each other within a simulation 
region. 

Overall, the model skill visualized in Figure 4 indic- 
ates that recorded variability both between and within 
European regions can be sufficiently well reproduced by 
the combination of migration and endogenous dynamics 
as formulated in GLUES. There is a significant correla- 
tion between the hindcastcd onsets and the reconstruc- 
ted onset in both the simulation with and without climate 
events (r^ = .37, r^ = .43, n = 39, p > .99, respectively). 
The simulation with climate events hindcasts the onset for 
most regions 162 years later than the average reconstruc- 
ted onset; the mean model bias is « 400 ± 1000 years; the 
difference to the simulation without climate events (me- 
dian onset 282 years earlier than the reconstruction) is 
statistically not significant. 

On average, climate events delay the onset by 461 ± 
179 years. Small delays (minimum 65 years) occur in re- 
gions where climate events are temporally separated from 
the transition, while long delays (maximum 1150 years) oc- 
cur in regions where several climate events occur shortly 
before and in the initial phase of the transition. Climate 
events during the transition tend to prolong the duration 
of the transition by 50-100 years but this trend is not stat- 
istically significant. 

4. Discussion 

4.I. Stagnation lines in model and archaeological evidence 
The Global Land Use and technological Evolution Sim- 
ulator is able to hindcast a realistic spatiotemporal pattern 
of the introduction of farming and herding into Europe 
between 7000 and 3500 simBC. Simulated transitions to- 
wards agropastoralism compare well to a large dataset of 
radiocarbon dated Neolithic sites. In the simulation, as 
well as in the data, agropastoralism did not expand uni- 
formly, but rather in periods of rapid spread interrupted 
by periods of spatial stagnation — but local intensification. 
This rapid Neolithization, for example, applies to the ex- 
pansion from Greece to the central Balkan in the 67th 
century simBC, which is followed by a 400a period of re- 
lative stagnation. A very similar pattern is hindcastcd for 
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Figure 3: Timing of the transition to agropastoralism in Western Eurasia. The simulated transition (background pastel shading) is contrasted 
with the radiocarbon ages of Neolithic sites from Pinhasi et al. (2005, solid color triangles). Bold lines indicate regional stagnation periods 
where the onset lag between neighboring regions is at least 500 years. 
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Figure 4: Reconstructed timing of agropastoralism versus simulated 
onset using a constant or a fluctuating climate forcing (grey and black 
circles, respectively). The lower left cluster of outliers represents 
Fertile Crescent and Anatolian founder regions. 



the Linearbandkeramik-like Neolithization in the 55th and 
54th century calBC, and for the relative stagnation before 
the onset of a Funnelbeaker-like cultural complex further 
north, in the 49th to 46th century. Other stagnation lines 
identified in the model occur around the Levant and along 



the Pyrenees, and divide the southern from the northern 
Iberian peninsula. 

These stagnation periods were archaeologically recog- 
nized in Guilaine's (2003) 'hypothese arythmique'. Based 
on this hypothesis, Rassc (2008) identified stagnation lines 
by analyzing isolines in the European Neolithic radiocar- 
bon record: he also found long duration stagnation lines 
running through Iberia and the Pyrenees, and one cross- 
ing Anatolia; shorter stagnation lines were identified in 
Greece, one bisecting the Balkan, another running along 
the Carpathian mountains, and one along the Alps. Ex- 
tending the earlier work, Schier (2009) emphasized a major 
stagnation line along the northern edge of the European 
loess belt, separating the Linearbandkeramik and Funnel- 
beaker cultures between 5100 and 4400 calBC. 

Dispersal models of the European Neolithic create spa- 
tial structure in the simulated Neolithization pattern from 
geographic and topographic constraints (e.g. Davison et al., 
2006; Ackland et al., 2007; Galeta et al., 2011). We add 
to this infrequent climate excursions emerging as sudden 
decreases of the natural productivity and diffusion rates 
of traits and people. This diffusion not only depends on 
the topography and geography, but again on (simulated) 
people and their technology. Compared to those dispersal 
studies, our model draws a more heterogeneous picture 
of the regional transitions to agriculture; it is the only 
model which is capable of representing stagnations in the 
expansion. We are thus one step further in establishing a 
simulation of the deterministic pattern of the Neolithic in 
Europe. Where this deterministic pattern fails to repro- 
duce archaeological information, we can in the future more 
precisely identify times and regions where human agency 



played a dominant role. 

Before the discrepancies can be exploited scientifically, 
several caveats about the simulation have to be considered: 
(1) the model deviation is larger than the uncertainty as- 
sociated with radiocarbon dates for most sites. (2) each 
simulation region covers a large area which cannot be fully 
represented by individual sites, such that the scale dif- 
ference introduces an additional comparison uncertainty; 
(3) there is a a spatial bias in the data: the number of 
sites per simulation region varies, and there is good cov- 
erage along the transect from the Levant to northwestern 
Europe, but few or no information on eastern Europe, on 
the Iberian peninsula, and in North Africa. 

4-2. Role of climate events 

Reconstructed climate events episodically depressed pop- 
ulation density and to a lesser extent technological stages. 
The impact of events on the timing of the Neolithic sim- 
ulation is largest whenever they occurred before or during 
the first stage (between 10-50% agropastoralists) of a local 
transition: In 29 regions (out of the 53 simulation regions 
in Europe), a climate event occurred during the transition. 
On average, the transition was delayed by 50-100 years 
(w.r.t. a simulation without events). In three cases, a re- 
gionally emerging agropastoral life style reverted back to 
hunting and gathering during a climate event. 

Not everywhere, however, where a climate event coin- 
cides with the transition, there is a delay; and vice versa, 
there are delays in regions which did not experience cli- 
mate events during the transition to agriculture, such as 
in most parts of western central Europe. All individual re- 
gions are closely connected by neighbor and remote trade 
relations: impacts, such as those of climate events, may 
have had far-reaching consequences beyond the locally af- 
fected region (Weninger et al., 2005). 

As evident from comparing this study with the one 
by Lemmen et al. (2011), climate fluctuations hindered 
the development in Greece and the Balkan (which is now 
more realistic) , while the temporal gaps between the cent- 
ral Balkan, the Linearbandkeramik culture, and the Fun- 
nelbeaker regions become more pronounced, which is also 
more consistent with archaeological evidence (Pinhasi et al., 
2005; Barker, 2006; Schier, 2009). Compared to the field 
evidence, the simulation shown here is too late in repro- 
ducing the transition to farming in France. Though mi- 
gration is not enhanced by climatic triggers in continental 
Europe, it docs increase so for Great Britain, where the 
transition occurs as early as France thanks to increased 
immigration. 

The differences we obtain between simulated timing 
and the radiocarbon age of related sites (Figures 3) in- 
dicate a partial improvement over the simulation without 
climate events shown by Lemmen et al. (2011): the overall 
bias of the simulation increased slightly (the event for- 
cing led to an overall delay of the onset) , but the variab- 
ility in the data could be redrawn more precisely. Des- 
pite the non-significant regressions between reconstructed 



and hindcasted onsets and durations of the Neolithic, the 
scatter alignment is close to the 1:1 line — apart from a 
bias in the onset hindcasts for the Levante regions (Fig- 
ure 4); neither regionally clumped chronologies nor model 
assumptions can be expected to be quantitatively precise. 

Overall, uncertainties in the simulated timing, the ra- 
diocarbon dates, and in the regional up-scaling are yet 
so large, that the differences between the simulation with 
and without climate events are not statistically signific- 
ant for Europe. Yet, the reconstruction of climate events 
improved the hindcasted Neolithic transition, both in tim- 
ing and regional transition durations. On the one hand, 
this points to an effective albeit limited sensitivity of soci- 
etal dynamics to environmental shifts. On the other hand, 
this result supports the notion of a large resilience of cul- 
tures, just like Grosman (2003) found for the Late Natufi- 
ans. Agency, the capacity of cultural complexes to adapt 
and preserve their life-styles still appears as the dominant 
model for understanding past human ecodynamics. 

This short paper is only a start in the direction of tack- 
ling the climate-culture link on large scales, and it gives 
a perspective of what could be achieved by bringing pa- 
laeoclimate data, cultural modeling, and chronologies to- 
gether. The idealization of climate events we present here 
needs to be discussed and evaluated further, likewise the 
(few) assumptions made by the model. In this attempt, 
the increasing availability of local chronologies will stimu- 
late regional-scale simulation studies for Europe and other 
continents. 

5. Conclusion 

We presented a spatially explicit mathematical model 
of the Neolithization of western Eurasia from 7000 BC to 
3500 BC. Our model incorporates endogenous sociotech- 
nological dynamics, as represented by the adaptation of 
characteristic population traits and their interaction with 
demographics and changing environments. The simula- 
tion improved in hindcasting the Neolithic expansion after 
integrating quasi-realistic climate events, especially with 
regard to stagnation phases visible also in field data. This 
outcome helps to confine the relevance of climate variabil- 
ity in infiuencing past human ecodynamics and suggests a 
dominance of endogenous over exogenous control factors. 
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Supporting material 

In the supporting online material, a movie of the ex- 
pansion of agriculture (Movie SI), the full table of glob- 
ally distributed palaeoclimate time series (Table SI), and 
a more detailed description of GLUES are available. The 
simulated data have been permanently archived on and 
are freely accessible from PANGAEA (Data Publisher for 
Earth & Environmental Science) as a netCDF dataset with 
reference doi:10.1594/PANGAEA.779660. GLUES is free 
and open source software and can be downloaded from 
http://glucs.sourceforge.net. 
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